Scaling of Charged Particle Production in d+Au Collisions at y/s 
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The measured pseudorapidity distributions of primary charged particles over a wide pseudo- 
rapidity range of \r\\ < 5.4 and integrated charged particle multiplicities in d + Au collisions at 
^/s NN = 200 GeV are presented as a function of collision centrality. The longitudinal features of 



-Au collisions at 



200 GeV are found to be very similar to those seen in p + A collisions at 



lower energies. The total multiplicity of charged particles is found to scale with the total number of 
participants according to Nj?j^ u = ^(N part }N^ and the energy dependence of the density of charged 
particles produced in the fragmentation region exhibits extensive longitudinal scaling. 

PACS numbers: 25.75-q,25.75.Dw,25.75.Gz 



Charged particle multiplicities have been studied ex- 
tensively in high energy collisions because of the intrinsic 
interest in understanding the production mechanism. 
More recent interest comes in the context of searching 
for and studying new forms of matter that are expected 
to be created in heavy ion collisions at ultrarclativistic 
energies. A key quantity that contains information about 
the longitudinal aspects of the multiparticle production 
process, and that has provided valuable input for discrim- 
inating between phenomenological models in the past 
0, is the rapidity distribution of identified particles. 
When particle identification is not available, the almost 
equivalent pseudorapidity distribution of charged parti- 
cles suffices. For this reason, such distributions have 
been studied in detail in hadron + proton fj, hadron 
+ nucleus and nucleus + nucleus Q collisions. 

Since the first collisions were achieved at the Rela- 
tivistic Heavy Ion Collider (RHIC), the PHOBOS col- 
laboration has obtained extensive information on pseu- 
dorapidity distributions of charged particles produced 
in Au + Au collisions for nucleon-nucleon center-of-mass 
energies (y / s^~Tv) between 19.6 and 200 GeV and over a 
large range of collision geometries. We have observed, 
for example, that the particle density in the midrapidity 
region changes smoothly as a function of ^/s NJV [5j and 
that the total charged particle production scales linearly 
with the number of participants 6] . Further detailed ob- 
servations of the shape of the pseudorapidity distribution 
show a scaling according to the "limiting fragmentation" 



hypothesis 0. The study of a simpler system such as 
d + Au at the same energy as Au + Au is essential 
to gain insight into which aspects of the data follow 
from the initial interacting states or general properties 
of the hadronic production process, and which are the 
consequence of the very different environments created 
in d + Au and Au + Au collisions. 

In this Letter we present the results of detailed mea- 
surements of the pseudorapidity distributions of primary 
charged particles, dN c h/d?y, as a function of collision 
centrality in d + Au collisions at y/s NN = 200 GeV Q 
over a wide pseudorapidity range of \r)\ < 5.4. The 
pseudorapidity, ?y, is defined as rj = — ln[tan(0/2)], where 
8 is the emission angle relative to the direction of the 
deuteron beam. The results for d + Au collisions are 
compared to Au + Au collisions and inelastic p(p) + p 
collisions at t/s^ = 200 GeV as well as to p(7r + , K + ) 
+ nucleus collisions at lower energies. 

The data were obtained with the multiplicity array of 
the PHOBOS detector 9] at RHIC. The array consists 
of an octagonal barrel of silicon detectors, the Octagon, 
surrounding the interaction region in an approximately 
cylindrical geometry covering \rj\ < 3.2. This array is 
augmented by two sets of three annular silicon counter 
arrays, the Rings, along the beam pipe far forward and 
backward of the interaction region (3.0 < \rj\ < 5.4). This 
array is identical to that used in our study of Au + Au 
collisions 0. The setup also includes two sets of 16 
scintillator counters (3.0 < \rj\ < 4.5) which were used 
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FIG. 1: Measured pseudorapidity distributions of charged 
particles from d + Au collisions at JIT^ 



200 GeV as a 



function of collision centrality. Shaded bands represent 90% 
confidence level systematic errors and the statistical errors are 
negligible. The minimum-bias distribution is shown as open 
diamonds [Tof. 



in the primary event trigger and in the offline event 
selection. 

The centrality determination was based on the ob- 
served total energy deposited in the Ring counters, ERi ng , 
which is proportional to the number of charged particles 
hitting these detectors. The choice of this centrality 
measure was based on extensive studies utilizing both 
data and Monte Carlo (MC) simulations that sought to 
minimize effects of auto-correlations on the final dN c h/d?7 
result. These effects can be significant when using other 
centrality measures |l0j , resulting in enhancements (sup- 
pressions) in the reconstructed midrapidity yields of up 
to ~ 30% for central (peripheral) collisions. The MC 
simulations used in the study included both HIJING 
[rll | and AMPT ^3] event generators coupled to a full 
GEANT P3 simulation of the PHOBOS detector. 

Four additional centrality measures, discussed in Ref. 
[To|. were created in order to study the detailed effects 
of auto-correlation biases. Ratios of the reconstructed 
dN c h/d?7 distributions obtained from the five centrality 
measures for data and, independently, for the MC simula- 
tions were found to agree, giving confidence in the entire 
methodology. This information, together with knowledge 
of the unbiased MC simulated "truth" distributions, pro- 
vided a clear choice of the centrality measure based on 
the Ring detectors as that which yielded the least bias 
on the measurement. It is important to note that this 
study only provided guidance with respect to the choice 
of Efjjng for the experimental centrality measure, and the 
final experimental dN c h/dry results do not rely in any way 



on the detailed shape of the dN c h/d?y distributions from 
the MC simulations. 

The multiplicity signals of Efu ng were divided into five 
centrality bins, where each bin contained 20% of the 
total cross section. For this to be done correctly, the 
trigger and vertexing efficiency had to be determined for 
each bin. Knowledge of the efficiency as a function of 
multiplicity allowed for the correct centrality bin deter- 
mination in data as well as the extraction of the corre- 
sponding efficiency-averaged number of participants. A 
comparison of the data and the MC simulations yielded 
an overall efficiency of ~ 83%. 

Results of the Glauber calculations implemented in the 
MC were used to estimate the average total number of 
nucleon participants, (Np art ), the number of participants 
in the incident gold, (Np art ), and the deuteron, (Np art ), 
nuclei, as well as the number of binary collisions, (N co ii), 
for each centrality bin (see Table QJ. 

The details of the analysis leading to the measurements 
of dNch/dr; can be found in Ref. [l4|. The measured 
dN c h/d?7 was corrected for particles which were absorbed 
or produced in the surrounding material and for feed- 
down products from weak decays of neutral strange par- 
ticles. Uncertainties in dN c h/dry associated with these 
corrections range from 6% in the Octagon up to 28% in 
the Rings. These uncertainties dominate the systematic 
errors. 

Figure ^ shows the pseudorapidity distributions 
of primary charged particles for d + Au collisions 
at y/s NN = 200 GeV in five centrality bins and for 
minimum-bias events. A detailed discussion of our 
minimum-bias distribution can be found in Ref. E3. As 
a function of collision centrality, the integrated charged 
particle multiplicity in the measured region (|?j| < 5.4) 
and the estimated total charged particle multiplicity 
extrapolated to the unmeasured region using guidance 
from the shifted p+nucleus data (see Fig. |2J) are pre- 
sented in Table |I| The centrality bins 0-20% and 80- 
100% correspond to the most central and the most pe- 
ripheral collisions, respectively. The pseudorapidity is 
measured in the nucleon-nucleon center-of-mass frame; a 
negative pseudorapidity corresponds to the gold nucleus 
direction. For the most central collisions, the mean r\ 
of the distribution is found to be negative, reflecting 
the net longitudinal momentum of the participants in 
the laboratory (NN) frame. For more peripheral col- 
lisions, the mean rj tends to zero as the distribution 
becomes more symmetric. For measurements of d+Au in 
the nucleon-nucleon center-of-mass system the Jacobian 
between dN c h/dy and dN c h/d?7 naturally produces the 
"double- hump" structure in dN c h/d?7 even if there is no 
structure in dN c h/dy. 

Now, we compare our d + Au results with p + A data 
obtained at lower energy, and discuss the energy and 
centrality dependence of the data. Figure |3 compares 
dN c h/d77 distributions of d + Au to p + Emulsion (Em) 
collisions at five energies 0, in the effective rest 
frame of both the projectile "beam" (a) and target (b). 
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FIG. 2: Panel a): Comparison of dNch/dr; distributions for 
d + Au collisions at ^/s NN — 200 GeV to p + Em collisions 
(sum of shower and gray tracks) at five energies mm. The 
rj measured in the center-of-mass system has been shifted to 
f] — ybcam in order to study the fragmentation regions in the 
deuteron/proton rest frame. Panel b): similar to panel a) 
but shifted to r\ + ytargot in order to study the fragmentation 
regions in the gold/Emulsion rest frame. 



For p + Em data the pseudo-rapidity distributions rep- 
resent the sum of shower and gray tracks. It should be 
noted that rj is measured in a different reference frame 
for d + Au and p + Em. This means that, compared to 
d+Au collisions, the p + Em pseudorapidity distributions 
are suppressed by the Jacobian for rj + ytarget ~ 0. The 
50-70% centrality bin of d + Au collisions was selected 
in order to match as well as possible Np art /Np art to 
N p^t/ N part where (N part (pEm))=3.4. The relative nor- 
malization of the dNch/dry for d + Au and p + Em colli- 
sions requires a ratio of (N par t(dAu))/(N par t(pEm)) = 1.6, 
such that the data correspond to the same number of 



30 
25 
20 
15 
10 
5 




a) 



200 GeV d + 
200 GeV d + 
200 GeV Au 
200 GeV p + 
200 GeV p + 



Au 

Au Min-Bias 
+ Au 

p inelastic 
pNSD 



10 



10' 



< Npart ) 



II 
< 
DC 



200 GeV d + Au 



O 200 GeV d + Au Min-Bias 



J 



b) 



1 = 



2^ 



E178 

K + + Pb 
7i + + Pb 
p + C 
p + Cu 
p + Pb 



6 



8 10 
< Npart > 



12 



14 



16 



pp 

ch ' 



FIG. 3: Panel a): Total integrated charged particle multiplic- 
ity per participant pair in d + Au, Au + Au @, p(p) + p in 
elastic |19l ] and p(p) + p NSD 20] collisions at the same en 
ergy, ^/s^ = 200 GeV. Panel b): The ratio R A = N ch /N 
where NJ?^ is the total number of charged particles for inelastic 
p(p) + p collisions, as a function of the total number of 
participant nucleons (N par t) for different collision systems. 
The 7r + +Pb, p+C, p+Cu and p+Pb collisions are for ^s NN = 
19.4, 13.7 and 9.69 GeV,and the K++Pb for ^/s^ = 13.7 
and 9.69 GeV from Ref. |2J , respectively. For the d + Au and 
Au + Au data the error bars indicate the systematic errors 
(90% C.L.). The dashed line represents the linear relation, 
Ra = f<N part ). 



nucleons interacting with the nucleus. A remarkably 
good agreement (limiting fragmentation) is observed in 
the fragmentation regions of the deuteron (gold) between 
d + Au and p + Em collisions at different energies. 
Furthermore, the overlap between the fragmentation re- 
gions of the deuteron (gold) and proton (Em) extends to 
lower | r\ | with increasing collision energy. In Ref. |17| we 
reported a similar comparison but for more central d+Au 
and p+Pb collisions. This extensive longitudinal 
scaling, also seen earlier in p(p)+p Q, P + A [l8j . 
and A + A collisions, suggests that the dominance of 
the two fragmentation regions is a common longitudinal 
feature of all multiparticle production processes. 

Figure |2K shows the total integrated charged particle 
multiplicity per participant pair for d + Au, Au + Aula] 
collisions, p(p) + p inelastic an d P(p) + P NSD |2(j 
all at the same energy, y/s~^ = 200 GeV, as a function 
of the total number of participants. The results show that 
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FIG. 4: Centrality dependence of the dNch/dr? ratio of d + 
Au collisions relative to that for inelastic p(p) + p collisions 
[2^| at the same energy. The arrows represent the (Np art ) of 
the most central and the most peripheral collisions. Typical 
systematic errors are shown for selected points (90% C.L.). 
The inset figure shows the variation of the ratio as a function 
of the number of gold participants, (Npa rt ), for four r\ regions 
in the gold direction. The dashed lines represent a linear fit 
to the data. 



the total charged particle multiplicity scales linearly with 
(Np ar t) in both d + Au and Au + Au collisions. They also 
indicate that the transition between inelastic p(p) + p 
collisions and Au + Au collisions is not controlled simply 
by the number of participants, as even very central d + 
Au multiplicity per participant pair shows no signs of 
extrapolating to the Au+Au results. Not only do we find 
that the total charged particle production in d+Au scales 
linearly with (N part ), but also that the scaling is energy 
independent and is the same in all hadron + nucleus 
collisions This is evident from Figure |3Jd where the 
ratio Ra = N c h/N^ is plotted as a function of (N part ) 
for a large variety of systems and energies. Here N cri is 
the integrated total charged particle multiplicity for d+ 
A, p + A [2lJ and N££ is for p(p) + p. It is this scaling, 
observed for the first time by Busza et al. ^ , which led 
to the wounded nucleon model of Bialas et al. |23| . 
Finally, in Figure 0] we compare for all centralities, 



the ratio of dN cri /d77 distributions observed in d + Au 
collisions relative to that for inelastic p(p) + p |24[ at 
the same energy. The inset figure shows the ratio as 
function of the number of gold participants, (Np art ), for 
four 77 regions in the gold direction. The results are 
consistent with a picture in which the production of 
particles with rapidity near that of the incident deuteron 
(gold) is approximately proportional to the number of 
deuteron (gold) participants. These results are consistent 
with lower energy p + A data pll I2H and, with the 
quark-parton model of Brodsky et al. [26j . 

In summary, we find that the longitudinal features of 
d+Au collisions at t/s^ — 200 GeV, as reflected by the 
centrality dependence of the pseudorapidity distributions 
of charged particles, are very similar to those seen in 
p + A collisions at energies lower by more than an order 
of magnitude. In particular, we find that in d + Au 
collisions the total multiplicity of charged particles scales 
linearly with the total number of participants, that the 
transition between the multiplicity per participant in 
d + Au and Au + Au collisions is not controlled simply 
by the total number of participants, and that the energy 
dependence of the density of charged particles produced 
in the fragmentation regions exhibits extensive longitu- 
dinal scaling. These results impose strong constraints on 
models of multiparticle production. 
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TABLE I: Estimated number of nucleon participants in 
the incoming gold ((Np" rt )) and deuteron ((Np art )) nuclei 
((N P art}= (Npart)+ (Npart)) and the number of collisions 
((N co ii)) as a function of centrality in d + Au collisions. 
The integrated charged particle multiplicity in the measured 
region (\r/\ < 5.4) and the estimated total charged particle 
multiplicity extrapolated to the unmeasured region (see text) 
are listed. All errors are systematic. 



Nfi<5.4 Nft t 



Cent. (%) (N^) (Npart) (Ncoii) 

O^Fj 13.5 ± 1.0 2.0 ± 0.1 14.7 ± 0.9 157 ± 10 167+jj 

20-40 8.9 ± 0.7 1.9 ± 0.1 9.8 ± 0.7 109 ± 7 115 

40-60 5.4 ± 0.6 1.7 ± 0.2 5.9 ± 0.6 74 ± 5 77 



id 

8 
+7 
5 

46 ± 3 48±3 



60-80 2.9 ± 0.5 1.4 ± 0.2 3.1 ± 0.6 

80-100 1.6 ± 0.4 1.1 ± 0.2 1.7 ± 0.5 28 ± 3 29±| 

Min-Bias 6.6 ± 0.5 1.7 ± 0.1 7.1 ± 0.5 82 ± 6 87tg 

50-70 3.9 ± 0.6 1.6 ± 0.2 4.2 ± 0.6 59 ± 4 62^! 



[1] N.N. Nikolaev, Sov. J. Part. Nucl. 12, 63 (1981); D. Horn, 
Phys. Rep. C 4, 1 (1972) and references therin. 

[2] J. Whitmore, Phys. Rep. 27, 187 (1976); G.J. Alner et 
al, Z. Phys. C 33, 1 (1986) and references therein. 

[3] W. Busza, Acta Phys. Pol. B 8, 333 (1977); S. Fredriks- 
son et al., Phys. Rep. 144, 187 (1987). 

[4] HR. Schmidt et al, J. Phys. G. 19, 1705 (1993). 



[5] B. B. Back et al. Phys. Re v. Let t. 88, 022302 (2002). 

[6] B. B. Back et al., Preprint nucl-ex/0301017 

[7] B. B. Back et al, Phys. Rev. Lett. 91, 052303 (2003). 

[8] Technically, the d + Au collisions were run at 
^/s NN — 200.7 GeV, matching the beam rapidities 
for Au + Au collisions at ^/s NN = 200.0 GeV, see 



5 



http:/ /www. agsrhichome.bnl.gov/RHIC/Runs/ 
[9] B. B. Back et ai, Nucl. Instr. Meth. A 499, 603(2003). 
[10] B. B. Back et al, Phys. Rev. Lett. 93, 082301 (2004). 
[11] M. Gyulassy et ai, Comp. Phys. Comm. 83, 307 (1994). 
[12] Zi-wei Lin et al, Phy. Rev. C 68, 054904 (2003). 
[13] R. Brun et al, GEANT 3.21, CERN Lib., W5013, 1994. 
[14] B. B. Back et al, Phys. Rev. Lett, 87, 102303 (2001). 
[15] Compilation of world data in p-Em collisions, see S. 

Fredriksson et al, Phys. Rep., 144, 187 (1987); I. Ot- 

terlund et al, Nucl. Phys. B 142, 445 (1978). 
[16] A. Abduzhamilov et al, Phys. Rev. D 35, 3537 (1987). 
[17] R. Nouicer (for PHOBOS Collaboration), J. Phys. G: 

Nucl. Part. Phys. 30, S1133 (2004). 



[18] C. Halliwell et al, Phys. Rev. Lett. 39, 1499 (1977); I. 

Otterlund et al. , Nucl. Phys. B 142, 445 (1978). 
[19] G.J. Alner et al, Phys. ' Rep. 154, 247 (1987) and 

H. Heiselberg, Phys. Rep. 351, 161 (2001). 
[20] G.J. Alner et al, Phys. Lett. B 167, 476 (1986). 
[21] J.E. Elias et ai, Phys. Rev. D 22, 13 (1980). 
[22] W. Busza et al, Phys. Rev. Lett, 34, 836 (1975). 
[23] A. Bialas et al. Nucl. Phys. B 111, 461 (1976). 
[24] G.J. Alner et al, Z. Phys. C 33, 1 (1986). 
[25] C. DeMarzo et al. Phys. Rev. D 29, 2476 (1984). 
[26] S. J. Brodsky et al, Phys. Rev. Lett., 39, 1120 (1977). 



